Prepared for submission to JHEP 



Spherically Symmetric Solutions in Massive Gravity 
and Constraints from Galaxies 



Stefan Sjors" and Edvard MortselP 

'^Department of Physics & The Oskar Klein Centre for Cosmoparticle Physics, 
Stockholm University, AlbaNova University Centre, SE-106 91 Stockholm, Sweden 

E-mail: stefan.sjors@fysik.su.se, edvard@fysik.su.se 

Abstract: In this paper, analytical solutions describing static and spherically symmetric 
sources in the decoupling limit of massive gravity are derived. We analyze the model 
parameter range and specify when a Vainshtein mechanism is possible. Furthermore, we 
use gravitational lensing and velocity dispersion data from galaxies to put constraints on the 
mass scale of the graviton. The result for the inverse graviton mass scale = h/{cmg), 
in units of the Hubble radius rn = c/Hq, is of the order Xg/rn ^ 0.01 — 0.02 at 95% 
confidence level. 

Keywords: modified gravity, massive gravity, gravitational lensing, velocity dispersion 



Contents 

1 Introduction 1 

2 Background 2 

3 Formalism 3 

3.1 Massless and massive linear gravity 4 

4 Model 5 

4.1 Action and equations of motion 6 

4.2 Static and spherically symmetric ansatz 7 

4.3 Solutions for vr' 8 

5 Method 10 

5.1 Velocity dispersion 11 

6 Data 11 

7 Results 12 

8 Conclusions 12 
A Solving for vr' 14 



1 Introduction 

The discovery of the late-time cosmic acceleration [1, 2] has drastically changed our pic- 
ture of the universe. In the current understanding, the energy content of the universe is 
dominated by a small cosmological constant, making the cosmological constant problem 
more pressing than ever. One possible direction towards an explanation of the cosmic 
acceleration and perhaps also a route to a solution to the cosmological constant problem 
might be offered by modifications of general relativity (GR) at large distances. As of today, 
there are many ideas of large distance modifications of GR, some more prolific than others. 
One simple idea of large distance modifications, one which we pursue in this paper, is the 
possibility of adding a small mass nig to the graviton. Needless to say, there are many 
other interesting scenarios such as the brane world of Dvali-Gabadadze-Porrati (DGP) [3] 
and scalar-tensor theories such as Chameleon [4] and f{R) theories [5]. 

These theories have in common that not only do they modify gravity at large dis- 
tances, but they also introduce small modifications to gravity around local sources. These 
deviations must necessarily be kept minimal so that the theory passes the very stringent 
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solar system tests. Here, non-linear effects play a crucial role in suppressing the deviations 
from GR. For example, in massive gravity (MG) and the DGP model, the recovery of GR 
in regions of strong fields is realized by the so-called Vainshtein mechanism [6], while in 
Chameleon theories, the Chameleon mechanism recovers GR in regions of high density [4]. 
Even though the deviations from GR are kept minimal, all theories give specific predictions 
for the deviation depending on the non-linear realization of the recovery of GR. This gives 
us an opportunity to constrain theories of modified gravity using dynamical measurements. 

The deviation from flat space induced by a static spherically symmetric source is in gen- 
eral described by two functions; the gravitational potential $ and the spatial curvature ^. 
In GR it follows from Einstein's equations that these two functions are the same. Whether 
they are the same in theories of modified gravity is a question of the specific dynamics. 
Probing the two potentials using local measurements we can constrain deviations from GR. 
The two potentials are probed using different methods. Massive non-relativistic observers 
are sensitive to the gravitational potential whereas massless relativistic observers are 
sensitive to the spatial curvature through the combination + ^'). 

In this paper we use dynamics of galaxies to constrain deviations from GR. We probe 
$ using velocity dispersion measurements of stars and the combination ^(^ + ^) using the 
gravitational deflection of photons. Gravitational lensing in the context of massive gravity 
was recently analyzed in [7]. 

This paper is organized as follows. In §2 we give a brief introduction to massive gravity. 
In §3 we present the general formalism and the relevant observables, while in §4 we discuss 
the specific model of massive gravity we consider. In §5 we discuss our method and in 
§6 we analyze the galaxy data. We present the results in §7 and conclude in §8. We also 
include an appendix § A where we give the details of the solutions to the equations of motion 
derived in §4. 

2 Background 

One of the oldest and simplest ideas of modifications of gravity at large distances is to 
give the graviton a small mass. The first successful step in this direction was taken by 
Fierz and Pauli [8] already in 1939 when they wrote down a consistent free (linear) theory 
of massive spin-2 particles. Consistent generalizations to interacting (non-linear) theories 
turned out to be much harder and examples were always plagued with ghosts. This lead 
Boulware and Deser to conjecture in 1972 [9] that there were in fact no consistent non-linear 
completions of massive gravity. Lately, the subject has spurred renewed interest when a 
family of actions of non-linear massive gravity was written down in [10] that was shown to 
be ghost-free to quartic order in non-linearities, in the so-called decoupling limit [11, 12]. 
Suggestions how to complete the action has since then appeared in the literature [13] and 
it was hoped that these actions would be ghost-free to all orders. It was recently shown 
by Hassan and Rosen [14-16] that a certain two-parameter family of non-linear massive 
gravity actions were indeed ghost-free. Thus it is exciting to consider the phenomenology 
of massive gravity and specially the explicit implementation of the Vainshtein mechanism. 
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At the heart of the observational differences between massless and massive gravity hes 
the van Dam-Veltman-Zakharov (vDVZ) discontinuity [17, 18]. Naively, one would think 
that any effects of a graviton mass could only be observed at scales of the inverse mass, 
where the exponential Yukawa suppression would kick in. In fact, on top of this effect, 
modifications to the gravitational force law persist even at shorter distances. This lead 
vDVZ to predict a 3/4 difference in the angle of light deflected around the sun in MG 
compared to GR, also in the limit of zero graviton mass. Already at that time, such a big 
discrepancy rules out FP theory, as the note added in proof reads of [17]. 

If that was the end of the story, massive gravity would as of then be ruled out ob- 
servationally. But Vainshtein noticed in [6] that linear perturbation theory, that is FP 
theory, breaks down at a distance scale ry, the Vainshtein radius. The Vainshtein radius 
comprise an intermediate scale between the very small scale of the source, the gravitational 
radius rg = 2G^M / ^ and the very large inverse mass scale of the graviton, the reduced 
Compton wavelength of the graviton \g = h/{cmg). Vainshtein showed that outside the 
Vainshtein radius there exists a perturbation series in ry/r, with FP theory at first order. 
However, at the Vainshtein radius massive gravity becomes strongly coupled and at smaller 
scales there exists a small r/ry limit such that GR is recovered. Whether there existed 
a non-linear theory that interpolated in between the two regimes was an open question. 
Numerical work on the subject [19] shows that this indeed happens in certain scenarios. 

Intuition about the Vainshtein mechanism can be gained by going to the so-called 
decoupling limit of MG. The decoupling limit offers a calculable regime where we can 
explicitly interpolate between the regions outside and inside the Vainshtein radius, thus 
allowing for explicit predictions being made also at the intermediate scale of the Vainshtein 
radius. The decoupling limit comprise a double-scaling limit taking Mpi — t- oo and rUg — )• 0, 
with a certain scale A^ = Mpi rrig fixed. This limit is suitable for our applications since 
heuristically we can think of this limit as describing massive gravity at distances smaller 
than the reduced Compton wavelength of the graviton, Xg — )• oo, yet on scales large enough 
that a weak- field limit of general relativity is approximate, Rs — ?• 0. This should be a good 
approximation since after all gravity is not observed to be exponentially weak at any of 
the scales we observe today, i.e. we have to push the graviton mass to a very small scale 
and throughout this paper we work in a regime where all distance scales are much shorter 
than the reduced Compton wavelength of the graviton, that is r <C A^. Furthermore, we 
are always far outside the gravitational radius of any object r ^ rs where linear gravity 
is an excellent approximation. Numerical solutions to massive gravity in the decoupling 
limit has been considered in [20] and exact results away from this limit in [21-24]. 

3 Formalism 

The most general static spherically symmetric perturbation of flat space is governed by 
two functions^ 

ds^ = -(1 + 2$/c2)(dx°)2 + (1 - 2^/c2)dx2 , (3.1) 
^Assuming general covariance. 
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where the gravitational potential ^ and spatial curvature ^ are functions of the radial 
distance r = |x| alone, and as usual we put = (ct,x), ^ = 0, 1,2,3. In linearized GR 
the two potentials coincide $ = but in theories of modified gravity the potentials are 
no longer equal. Within the parametrized post-Newtonian (PPN) formalism, the deviation 
from Newtonian gravity is quantified in terms of the ratio of the two potentials For a 

constant ratio, this defines the parameter 7ppn = ^/^, with 7ppn = 1 in linearized GR. In 
general, deviations from GR are not captured by a single number but rather two functions. 

The two potentials can be probed using different methods. A massive non-relativistic 
observer is only sensitive to the gravitational potential $ while a massless relativistic 
observer is sensitive also to the curvature ^. Specifically, a massless observer is sensitive 
to the potential combination = |(<I> + ^'), known as the lensing potential. Thus 
complementary observations using massive and massless probes can access both potentials 
and allows us to quantify deviations from GR. 

More quantitatively, we analyze the above geometry, eq. (3.1), using a probe with 
energy- momentum tensor t^,^. For a point particle moving on a trajectory x = x(t), the 
energy-momentum tensor is given by 

r^'^(x')=c^5(3)(x'-x(t)). (3.2) 

The force felt by the probe is F(x) = — V[/(x) where f7(x) is the interaction energy 
between the probe and the background. For a probe that is minimally coupled to gravity 
the interaction energy is given by 

U{^) = -\j d^x' h^'^i^) r^,(x') , (3.3) 

where h^y = g^y — rj^i, is the deviation of the metric, in eq. (3.1), from flat space. For 
a non-relativistic, nearly stationary particle of mass m and four-momentum = (mc; 0) 
we pick up /loo = — 2^>/c^ in the above and get for the interaction energy U = — m<I>. This 
reproduces the standard force law of a non-relativistic observer 

F = -mV$ . (3.4) 

For a massless particle = the energy momentum tensor is traceless t^'^ = and 
we are only sensitive to the traceless part of the metric perturbation which contains the 
combination $+ = i($ + indeed /jj|[ace-free ^ j^^^ _ ^rj^^h = -^+/c'^ ■ diag(3, 1, 1, 1). 
Taking for example p'^ = {E/c;0,0, E /c) the interaction energy becomes U = —2E^+/c^ . 

3.1 Massless and massive linear gravity 

As an example, consider a point source with mass M. In the weak-field limit of GR, the two 
potentials in eq. (3.1) are the same and equals the Newtonian potential <I>n = —G-^M/r 
of the point source 

$|gr = ^|gr = -^. (3.5) 
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Thus, in GR 7 = ^'/<I> = 1. As noticed by vDVZ, the situation is quite different in FP 
theory of massive gravity where^ 

$ FP = ^e-"/^« , ^' FP = ^e-''^^ , 3.6 

Sr Sr 

see for example the nice review [25]. At long distances, gravity is exponentially weak 
being cut-off by the graviton mass whereas at shorter distances, the gravitational force 
is effectively A massive observer probing short distances therefore experiences a 

Newtonian potential $|fp = —G'^M/r, albeit with a modified Newton's constant G'^ = 
4Gn/3. The degeneracy with GR is lifted using massless observers who are sensitive also 
to the curvature perturbation ^ and Fierz-Pauli theory predicts 7|fp = (^/$)|fp = 1/2. 
More strikingly put, since the lensing potential is the same in massive gravity as in massless, 
$+|fp = + ^)|fp = *^N! the prediction in terms of the modified Newton's constant 
of the deflection angle of a light-ray passing the perimeter of the sun is a factor of 3/4 
different in massive gravity than in massless, regardless of the graviton mass 

«, = ^^ = ^xi^. (3.7) 

r0 4 tq 

We choose to quantify the deviations from the force laws in GR defining two functions 
ecj>(r), evi/(r), such that defines GR, 

V$(r) = [l + e<i,(r)]V$N(r), (3.8) 
V^{r) = [l-e^{r)]V^^{r). (3.9) 

Notice that e$, e$ effectively changes the Newton's constant we would infer in a local mea- 
surement. For example, a non-relativistic observer would effectively experience a slightly 
varying Newton's constant G™(r) = [1 +e$(r)]GN while a relativistic observer would see 
a Newton's constant Gf^(r) = (1 + [£$(?") — e^j{r)]/2)G]^. 

Notice that in a phenomenologically viable theory of modified gravity, e$(r) and e^{r) 
necessarily vanish at short distances such that solar system constraints are evaded. For 
example, within the decoupling limit of massive gravity, e$(r) and e>i'(r) are slowly varying 
functions that vary from zero deep inside the Vainshtein radius, reaching one third at scales 
far outside the Vainshtein radius, i.e. 

lim e$(r), e4,(r) = , lim e$(r), ^^^(r) = 1/3 . (3.10) 

r/rv— >0 r/ry-^oo 

4 Model 

An intermediate step towards a ghost-free theory of massive gravity was taken in [10]. Here 
massive gravity was analyzed in the so-called decoupling limit and a certain two-parameter 
family of actions were shown to be ghost-free around fiat space. The decoupling limit 
comprise a calculable approximation of massive gravity and corresponds to taking the limit 



^Up to gauge equivalent terms. 
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Mpi — )• oo and nig — ?■ 0, with a certain scale = Mpi rrig fixed. This hmit is suitable for 
our applications since heuristically we can think of this limit as describing massive gravity 
at distances smaller than the reduced Compton wavelength of the graviton, — >■ oo, yet 
on scales large enough that a weak-field limit of general relativity is approximate, — )• 0. 

4.1 Action and equations of motion 

We now discuss the action of massive gravity in the decoupling limit. At the linearized level, 
the theory describes a massive spin-2 field with five degrees of freedom. In the decoupling 
limit, where the mass is taken to zero, the two helicity-1 modes become longitudinal (in 
the four-vector sense) and decouple from any conserved source. What remains is the action 
of the two helicity-2 modes and the one helicity-0 mode. Following [26] the action is given 
by^ 

^ = -^-H^-^S^^^H^fs + H^" + f-^xj^J + ^xj^'^ + ^H^'T,. . (4.1) 

The first term in eq. (4.1) comes from evaluating the Einstein-Hilbert action Mpj-^— gi? 
to quadratic order in the metric fluctuation M^^H^y = g^^ — rj^i,. Explicitly, the so-called 
Einstein operator S^l^ takes the form 

S^I^Hap = -\ {d^H^u - 2d^^dPH,)p + d^d^H - Tj^yd'^H + r^^.d^dpH'^'^) . (4.2) 

The second term in eq. (4.1) is a linear coupling of H^y and tt consistent with the Galilean 
symmetry vr — )• ir+c+h-x, forcing vr to enter in the combination 11^,^ = d^dyT^, together with 
gauge invariance H^^y — )• H^y — 29(^^,^), forcing X^i'^'^^ to be divergence free. Explicitly 

4i^=e/'"^^/p.n„^, (4.3) 

4'^ = e/^'^e/\n„^n,, , (4.4) 

x^} = e/'"^e/^'n,^n,^n.5 • (4.5) 

The term a H>"'xl]} gives a kinetic mixing between Hfj,^ and vr, and positivity of the kinetic 
energy gives a < 0. Apart from the sign, a is just a normalization of the field vr. The 
parameters (3, 7 are free model parameters. Setting a, /3, 7 = recovers weak-field GR. The 
third term in eq. (4.1) is the interaction energy density of source and gravity, as discussed 
in connection to eq. (3.3). 

From the above action, eq. (4.1), we can derive the metric equations of motion 

_ ^(1) , + + ^ T (A 6) 



and the tt equation of motion 

^Throughout §4, and only this section, we use natural units and put fi — c = 1. 



(4.7) 
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4.2 Static and spherically symmetric ansatz 

In this section we analyze the equations for a static, spherically symmetric, and pressure- 
free source described by a mass density p = p{r). The most general ansatz consistent with 
these symmetries can be put on the form 

ds^ = -[1 + 2$(r)] dt^ + [1 - 2^{r)] dx^ , (4.8) 
vr = 7r(r) , (4.9) 
V = diag(p(r), 0,0,0). (4.10) 

We now substitute the above ansatz into the equations of motion. Of the ten equations 
of motion of the metric, only the time-time and radial-radial component equations give 
independent equations. Together with the vr equation of motion we have three independent 
equations for three unknowns and the system should admit a solution. 

The two metric equations of motion read (using ' = to denote a radial derivative) 

= 2i|:; + " + i + F ^A-'f)' ' (4-11) 

^Pl'(^'-^') = -«'^'-|,^(v^')^ (4.12) 
while the vr equation of motion reads 

:0. (4.13) 



ar^{2^>' - + ^rvr'(^' - - ^^'{i^'f 



First, note that eq. (4.11) is a total derivative and can be integrated with the result 
*' = ^^ + ^vr' + ^i^ + ^i^ (414) 



16vrM2j r2 Mpi MpjAS r MpjAS 



where M{r) denotes the integrated mass inside the sphere of radius r 

M{r) = 4.TT [ dRR^p{R). (4.15) 
Jo 

This result is a manifestation of 'Gauss' law' or the 'Shell theorem': The force experienced 
at a radius r only depend on the integrated mass inside that radius, and the force is 
equivalent to that of a point particle of mass M situated at the origin. This property is 
an artifact of decoupling limit of MG, where the exponential Yukawa decay is pushed to 
infinity. Certainly, this is only a property of l/r^ forces and in the full theory this result 
does not hold. 

We now discuss briefly the two sub-cases GR: a,/3,7 = and FP: a < 0, /3,7 = 
before going on to the general case. 

GR case: From eq. (4.12) we see that in the GR limit, putting a,/3, 7 = 0, then indeed 
^' — ^i' = 0. Using eq. (4.14) we reproduce the standard Newtonian force law 
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where we in the last step identified the Newtonian force <I>^ = G^M{r) /r'^ , with Newton's 
constant^ Gn = (IGvrMpi)"^. 

FP case: In the case a < 0, /3,7 = 0, integrating eq. (4.13) forces 2^*' — <!>' = 0, (setting 
the integration constant to zero). Then using eqs. (4.12) and (4.14) we find 

3 2 (y. 1 

which reproduces eq. (3.6), (which also verifies the correct choice of integration constant). 

General case: Away from the GR limit we parametrize the deviations in the force laws 
using the two functions, e$ = ^*V*^n ~ ^ ^^"^ e<i- = 1 — ^''/^n; defined in eqs. (3.8) and 
(3.9). Using eqs. (4.12) and (4.14) we find 

a tt' 7 (tt')^ 

* Mpi MpiA3 r<^'^ MpiA6 r"^^'^ ' ^ ' 

Thus, the deviations are completely determined by the behavior of the vr field, or more 
specifically by powers of tt'. Substituting the expressions for and into eq. (4.13) for 
vr, we find a closed algebraic equation for vr'. Solving this equation for vr' then determines 
<!>' and ^' completely. Indeed, integrating eq. (4.13) (setting the integration constant to 
zero) and using the expression for ^' and $' we find a closed algebraic equation for vr' 

,.2^2 



"^ ^N + ^vr'+^^-^j (vr')^ 

^^±4a7 3 3r 

A^Mpi ' Ai2Mpir 



+ '^M^-X^(^?-0. (4.20) 



4.3 Solutions for vr' 



We now discuss the general features of the solutions to eq. (4.20), relegating the details to 
§A in the appendix. First, note that eq. (4.20) is a quintic polynomial equation in vr' and 
will therefore in general not allow for closed- form solutions^. However, as we show in §A 
the quintic term is negligible in almost all of the parameter space [/3,7] and we can drop 
the (vr')^ term from eq. (4.20), thus obtaining a cubic equation which we solve analytically. 
For 7 = 0, we find two solutions depending on the sign of 

B = ^lc?. (4.21) 

For sign(i?) = +1, the solution blows up at infinity and does not satisfy the right boundary 
conditions. For sign(i?) = —1, we find a solution with all the right boundary conditions 
as specified in eq. (3.10), allowing for a Vainshtein mechanism. The analytic solution for 



''Note our normalization of Afpi which is different from other normahzations often used in the literature. 
^Note that as a polynomial in r, eq. (4.20) is actually a quartic and r — r{Tv') can be obtained in a 



closed form. 
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B < 0, C = 



B = 0, C > 




Figure 1. Left: Results for C = and B < 0. Right: Results for B ^0 and C > 0. We plot the 
deviations in the force laws e$ = ^' /^'j^ — 1 and = 1 — 'i'' /^'^, together with deviation in the 
lensing potential $+/$n = 1 + ""^ 2"'" ^'^'^ ^ 1 = ^ 1- "^^"^ solutions exhibit the 

Vainshtein mechanism, eq. (3.10), i.e. lim^/r^^-i-o £*('')) ^*(^) = 0) lirnr/rv-i-oo £*('') = 1/3. 

Furthermore, the lensing potential is essentially that of GR everywhere with small corrections at 
r ~ ry, i.e. 1 + ~ 1 almost everywhere. 



7 = and with sign(S) = —1 is presented in eq. (A. 8) and plotted in Figure 1. The 
figure clearly demonstrates the Vainshtein mechanism; interpolating between GR at short 
distances r <C ry and FP at long distances r » ry. 

For /3 = we likewise find two solutions depending this time on the sign of 

C = -f/a^. (4.22) 

For sign(C) = —1, the solution blows up while for sign(C) = +1, there is a solution with all 
the right boundary conditions, as specified in eq. (3.10). The analytic solution for (3 = 
with sign(C) = +1 is presented in eq. (A. 11) and plotted in Figure 1. 

The solution for non-zero /3,7 is presented in eqs. (A. 13) and (A. 20). Eq. (A. 13) 
covers the range C > 0, B < while eq. (A. 20) covers the range C > 0, < i? < -Bmax- 
At -Bmax = ^5 — VT3V2C the quintic term gives important contributions and a singular 
behavior is developed beyond -Bmax, see Figure 3 in §A. We conclude that the quintic 
equation (4.20) allows for solutions with a Vainshtein mechanism, satisfying eq. (3.10), 
in the parameter range C > with B < -Bmax- Furthermore, the solutions to the cubic 
equation, neglecting the quintic term, approximates the solutions to the quintic equation 
within 0(10"^) in all of the parameter space C > 0, B < i?max, except at ~ -Bmax where 
the quintic becomes important, and there are up to 40% corrections to the cubic solution. 
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5 Method 



We now turn to the observational consequences of the above considered models. As 
explained in §3 we use the fact that a non-relativistic and a relativistic observer expe- 
rience different effective gravitational couplings, G™(r) = [1 + e<i.(r)]GN and Gj^{r) = 
(1 + ^[e^{r) — e^r(r)])GN respectively. Thus, deviations from GR can be probed by com- 
paring the masses of galaxies as measured from gravitational lensing and the velocity 
dispersion of stars within the galaxies. In practice, we compare the observed velocity 
dispersion within galaxies with the theoretically expected velocity dispersion assuming a 
model for the luminosity and mass distribution where the normalization of the mass is 
set by the image separation of the lensed images. The theoretically inferred velocity dis- 
persion will thus include the modified gravitational dynamics for both non-relativistic and 
relativistic observers. 

In the following, we will to a large extent follow the methods outlined in (author?) 
[27] (see also (author?) [28]). We will assume that the mass densities p and the luminosity 
densities u of the lensing galaxies can be written as power laws 

p = por~" , V = vqt'^ . (5.1) 

The spherical mass inside radius r is then 

Mir) = 47r T dRR^p(R) = i^^r^"" , (5.2) 
Jo 3 - a 

and the projected cylindrical mass traced by lensing within radius R 

27r 3/2A(a) po 

3 — a 

where A(x) is the ratio of gamma functions \{x) = r(^^)/r(|). Specifically, the projected 
mass within the Einstein radius Re 

M„„,(iJ,)=?!^!^M«fl|-.. (5.4) 

Since the angular separation of lensing images provides a good approximation to the cor- 
responding Einstein angle 9e = Re/ Di we can use lensing data to fix the normalization of 
the galaxy mass using 

"^-'<^-' = 4Gfk) iS^''- <^-^' 

where Di,Ds and Dig are angular diameter distances between the observer and lens, the 
observer and source, and the lens and source, respectively. Note in eq. (5.5), we use the 
appropriate effective Newton's constant G'^{r) being relevant for light. Furthermore, when 
deriving eq. (5.5) we have used that G^(r) is slowly varying such that when evaluating 
the light defiection line-of-sight integral, G]^{r) can be approximated by GJ^{Re)- 

We note that the ratio of the Einstein radius Re and the Vainshtein radius ry is given 

by 

3 /ReV DisDf . 
Pe= [ — ] = -fVvr^^E , (5.6) 

which is ^ 1 unless X„ <^ Di. 



rv J DgXl 
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5.1 Velocity dispersion 

Equations of stellar hydrodynamics give for the radial velocity dispersion 



a'Ar) = I G^\rHr)M{ryf^-' , (5.7) 

where (3 = 1 — (cJi/cJr.)^ is the (constant) velocity anisotropy of the system^ and Gf^{r) is the 
effective gravitational coupling as felt by massive observers. Since G^{r) = [1 + e$(r)]GN, 
we can write the velocity dispersion as a sum of the familiar Newtonian expression and a 
term depending on deviation in the force law e$(r) 

^r(0=4W+<(0, (5.8) 
where cr^{r) is the standard Newtonian result and o"^^(r) is given by 

poo 

< (0 = I e^{r).{r)M{ry^-^ . (5.9) 

The observed velocity dispersion, cr*, is effectively luminosity- weighted along the line 
of sight and over the effective spectrometer aperture. This averaging can be expressed as 



^ dR R wiR) JZ, dZ u{r) ( 1 - (3$) a^^r) 



= R ^ , (5.10) 

J,'^-dRRw{R)fZdZi^{r) 

where Z'^ = — R? and 

w{R) ^ e-^'/2a2,- , (5.11) 

is the aperture weighting function in which we will use a median seeing of CTatm = 1-4" and 
a cut-off radius of iimax = 3" [29]. In eq. (5.10), the factor ^1 — takes into account 

how the radial and tangential components of the velocity dispersion project along the line 
of sight. 



6 Data 



We use the data for the full sample of 131 strong gravitational lens candidates observed with 
the Advanced Camera for Surveys (ACS) aboard the Hubble Space Telescope (HST) by 
the Sloan Lens ACS (SLACS) Survey [29]. The lensing foreground galaxies are primarily of 
early-type morphology, with redshifts 0.05 < z < 0.5 and velocity dispersions 160km/s < 
Cobs < 400km/s. The lensed background galaxies have redshifts 0.2 < z < 1.2. The 
SLACS lens sample is statistically consistent with being drawn at random from a parent 
sample of comparable Sloan Digital Sky Survey (SDSS) galaxies, although our analysis 
does not depend on this property. Out of these systems, we use a subsample of 53 systems 
with elliptical lensing galaxies where singular isothermal ellipsoid gravitational lens models 

®Note that following conventions, we use the same notations a and j3 in describing the density profile 
and the velocity anisotropy as used as coefficients in the action eq. (4.f). 
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can successfully be fitted to the imaging data and we have reliable velocity dispersion 
measurements. 

For the present work, the key observables in each system are the redshifts and the 
stellar velocity dispersion of the lens galaxy (as measured from SDSS spectroscopy) and 
the Einstein angle of the strongly lensed image of the background galaxy (as measured from 
ACS imaging). We also make use of the effective radius measured from ACS imaging data 
to obtain individual estimates of the luminosity profile of the lensing galaxy. Following 
[29] , we adopt a 7 % velocity dispersion uncertainty and a 2 % Einstein radius uncertainty. 
We also add an additional 5 % to the velocity uncertainty given that we do not expect the 
observed velocity dispersion to perfectly match the one expected for a Singular Isothermal 
mass profile [30]. Note that the choice of luminosity profile may potentially alter the results, 
e.g. going from a power law to a Hernquist profile [31]. 

7 Results 

We are now in a position to generate constraints on the massive gravity model by comparing 
the observed velocity dispersions dobs with the theoretically calculated velocity dispersions 
for the given galaxy and massive gravity model a^:. From the analysis in the previous 
section, we have 

a, = a4a,^,5,9E;Xg,B,C), (7.1) 

where the semicolon separates the galaxy parameters and observables {a, (3, 6, Be) from the 
global massive gravity parameters (A^, B, C) which we are seeking to constrain. The slope 
and anisotropy of the lensing galaxies, a and /3, are being marginalized over on an individual 
basis for each galaxy assuming a prior knowledge of a = 2.00 it 0.08 and /3 = 0.13 it 0.13 
(68 % confidence level) [27] . We also make use of the effective radii measured from ACS 
imaging data to obtain individual estimates of the luminosity profile power-law index, 5, by 
comparing the observed total luminosity to the luminosity within half the effective radius. 

In each fit, we fix the value of B and C and constrain the value of \g in units of the 
Hubble radius rn = c/Hq. Results are shown in Figure 2. The cosmological limits are 
rather insensitive to the values of B and C and are on the order of Xg/rn ^ 0.01 — 0.02 at 
95 % confidence level. 

8 Conclusions 

We have analyzed large distance modifications of gravity in the specific context of mas- 
sive gravity. By deriving spherically symmetric solutions, we are able to constrain the 
parameters of the model using comparison to galaxy observations. The theory allows for 
a Vainshtein mechanism is some parts of the parameter space. In terms of the normalized 
couplings of the cubic and quartic interactions B = and C = the theory ex- 

hibits a Vainshtein mechanism for C > and B < i?max = \/5 - V^V2C. For C7 < 0, 
general relativity is not recovered at short distances, and worse for B > i^maxi the solution 
have pathologies becoming 'multiple valued' at B = -Bmax and 'singular' &t B = -v/SC*. 
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B = -1, C = 



B = 0,C =1 




0.02 0.04 0.06 0.08 0.1 0.02 0.04 0.06 0.08 0.1 

B = -l,C = l B = 1,C = 2 




0.02 0.04 0.06 0.08 0.1 0.02 0.04 0.06 0.08 0.1 



Figure 2. Cosmological constraints using strong lensing data with B = — 1,C = [upper left 
panel], B = 0,C — 1 [upper right panel], B — —2, C — 1 [lower left panel] and B — 1,C — 2 [lower 
right panel]. The limits are rather insensitive to the values of B and C and are on the order of 
>'g/rH > 0.01 - 0.02 at 95% confidence level. 



Analyzing gravitational lensing data and velocity dispersion measurements of galaxies 
we put limits on the inverse graviton mass scale Xg/rn ^ 0.01 — 0.02 at 95% confidence 
level in regions where we do have a Vainshtein mechanism. This pushes the graviton mass 
closer to the Hubble radius. An alternative approach would be to use clusters of galaxies 
acting as gravitational lenses instead of galaxies. In this case, the mass estimate using 
massive observers would be derived using X-ray data to observe the local pressure gradient 
inside the cluster to infer the gravitational potential. Although clusters would potentially 
allow us to probe larger radii than galaxies, we note that such an analysis rely on the 
assumption of cluster thermal equilibrium. 
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A Solving for n' 

Let us first formulate the problem of solving eq. (4.20) using more natural variables. Since 
we are ultimately interested in the deviations in the force laws, eqs. (4.18) and (4.19), let 
us consider the following dimensionless combination 

" ^' (A.l) 



Mpi 

In the following we will impose the boundary conditions e(r) — )• 0, when r — )• 0, and 
e(r) — )• 1/3, when r — )• oo, such that we reproduce GR close to the source and FP theory 
far away from the source. We will also find it convenient to take out the strong coupling 
scale rv = (r^A^^/^, where rs = 2GnM, from the radial variable defining a dimensionless 
variable^ 

r 

p^-. (A.2) 

rv 

Finally, it is more natural to incorporate powers of a to the couplings /? and 7, defining 

B = ^, C=^, (A.3) 

such that the rescaling 7r(r) — )• A~"'^7r(r) and a — )■ Xa, f3 — )• and 7 — t- A'^7 (which leaves 
the action in eq. (4.1) invariant) leaves B and C invariant. 
Using these variables, eq. (4.20) for it' takes the form 

while the deviations in the potentials, eqs. (4.18) and (4.19) take the form 

C o 

£^ = e- j^e , (A. 5) 

We now consider solutions to eq. (A.4), first for the two subcases B ^ 0, C = and B = 0, 
C 7^ 0, then attacking the general case B,C 7^ 0. 



^Not to be confused with the galactic mass density p. 
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Case: B ^ 0, C = For C = 0, the quintic eq. (A. 4) is reduced to a cubic in e, which 
allows for explicit solutions. Absorbing the coupling B into the radial variable, defining 
X = \B\~^/^p, eq. (A. 4) is reduced to 



1 - 3e + sign(B) 



0. 



(A.7) 



x^ 2x^ 

There are three solutions to this equation, one everywhere real and two somewhere complex 
solutions. For positive B, i.e. sign(i?) = +1, the real solution blows up at infinity and we 
throw this solution away. For negative B, i.e. sign(i3) = — 1 we find a real solution which 
stays finite at infinity. The real solution for B < and C = is given by 



£{x) 



u 



1/3 
B ' 



x) + — 



1/3/ N 



Ax-' 



\B\-^^^p, B<0, 



where ub{x) = — 16x^ + + 4x^\/— IGx^^ — 16x^ + 4. The behavior at small and large 
X can be inferred directly from the solution, eq. (A. 8), or more easily from eq. (A.7). For 
small X the cubic term dominates while for large x the linear term does, and we find 

1 



lim e{x) = 2^/^x^ , lim e 

x— >0 5>oo 3 



(A.9) 



Thus in the limit of small x we reproduce GR while in the limit of large x we reproduce 
FP theory. 



Case: i? = 0, C 7^ In the case B = Q, eq. (A. 4) takes the form 

3^2 ^3 3^5 

1 — 3e — sign(C)— 7 — sign(C)— + 



4y 



16y 



12 



0, 



(A.IO) 



where y = |C|~^/®/9. We now argue that in this case, the quintic term is of no importance. 
First, we demand that at infinity, e goes to 1/3. Then clearly e^/y^^ is highly suppressed 
compared to the other terms for large y. Second, for small y we demand that e goes to 0. 
If the quintic should give any contribution at all, then e oc y^^/^ at the origin, in which 
case the quadratic jy^ oc y24/5-6 i^Jows up and dominates the behavior at the origin. 
Thus the quintic term is of no importance neither at the origin or at infinity. We then solve 
the cubic equation, neglecting the quintic term, and simply observe that on the solution 
the quintic is negligible for all y, with maximal corrections 0(10"^) around y = 1. The 
cubic equation has one everywhere real solution and two somewhere complex solutions. For 
negative C, i.e. sign(C) = — 1, the solution does not go to zero at the origin, and we throw 
this solution away. For positive C, i.e. sign(C) = +1 we find a well-behaved solution 



e{y) = \ 



u 



1/3 

■c 



iy) + 



1 - 



y = \c\ 



-1/6, 



c > 



(A.ll] 



where is now easy to see, either directly from 

the solution, eq. (A. 11), or from eq. (A.IO), that we indeed have 



2 o 

linie(y) = —y 
y^o V3 



lim e{y) = - 

y-^co 3 



(A.12) 
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Case: B,C ^ In the case of non-zero B and C the situation is more subtle. Let us 
first analyze the parameter range of C. Now, observe that for C 7^ 0, the small p behavior 
is dominated by the terms 1 and — in eq. (A. 4). Thus, we conclude that e 



/3C 



P 



when p — )■ 0, which implies C > 0. We can than analyze the parameter range of B given 
some C > 0. We split up the case in S < and i? > 0. 

Subcase: C > 0, B < For C > and B < the situation is very much the same 
as in the already considered above cases. We find a solution to eq. (A. 4) (neglecting the 
quintic term) with the right behavior at the origin and at infinity which simply interpolates 
between eqs. (A. 8) and (A. 11) 

1 



2(2C + 52) 



where we have defined 

»(P)=9(P) + '1(2C + B2)p»v7(rf, 



u 



'/Hp) 



v{p) = ABp^-C, 

w{p) = (-16C + 8S2)p6 _ 8BCp^ + C2 , 
f{p) = 16(4C - ^2)^12 + 8(2^3 - 9BC)p^ 

+ (45^ + 37C2 - 20B^C)p^ + UBC^p^ - , 
g{p) = 16{B^ - 6BC)p^ + 4{2B^ - B^C + 14^2)^6 ^ i2BC^p^ 



(A.13) 

(A.14) 
(A.15) 
(A.16) 
(A.17) 
(A.18) 
(A.19) 



Just as the case C 7^ 0, i? = we find that the quintic term contributes negligible on the 
solution and eq. (A.13) is indeed a good approximation to the full solution. 

Subcase: C > 0, B > We then try solving eq. (A. 4) for B > 0, neglecting the quintic 
term. This time, the solution to the cubic equation is only a good approximation for 

B < Smax, where S^ax = 



13v2C as we show below. For B < -Bmax we find a 
solution with the correct limiting behavior 

^'/'ip) + 7Mr^+vip) 



<P) 



2(2C+_B2) 
1 



2(2C+B2) 

where is the zero of /(p*) 
induces important corrections 



«l/3(p) 

gi7r/3^1/3(-p-) _ g-i7r/3 w(p) 



v{p) 
2 



0. 



ul/3(p) 

We now show that for B 



for p< p^ 
for p> p^ 



(A.20) 



-Bmax the quintic term 



Geometric formulation: Let us consider the problem in more geometric terms. The 
curve e = e(p) in the pe-plane is the zero locus of the following polynomial in p, e with 



moduli B, C 

p{e,p;B,C) 



1 - 3e + 3 



B 



C_ 
4^ 



B^ C 

+ ^ 



3C2 



(A.21) 



^p^ 4/9" J \2p^ ' P'^ J' ' 16pl2 

As we vary the moduli B, C the curve is changing shape and we might encounter singular- 
ities, i.e. points where the tangent vanishes, see Figure 3 for reference. To find potential 
singularities we analyze the system of equations 



p = 0, dpp = 0, deP = 0. 



(A.22) 
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B > 0, C = 1 




Figure 3. As we vary the model parameter B = while keeping C = fixed, the curve 
e — e{r/rv) develops a singular behavior at i? = VSC where the curve 'opens up', see pink dotted 
line above. Already before this happens the curve becomes 'multiple valued' at i? = \/5 — VT3\/2C, 
see the dashed curve which acquires three values for e{r/rv) in the region r/ry ~ 1.2. For B < 
\/5 — \/l3-\/2C' the curve is well-behaved, see orange solid curve. The above example highlights 
Vainshtein's original idea: Models with a well-behaved short and long-distances behavior might or 
might not analytically connect the two regions. 



Using standard methods we can find a Groebner basis where the variables p, e are elimi- 
nated, and the system in eq. (A. 22) is reduced to 

B'^ -15 B'^C + 360"^ = 0. (A.23) 

This equation has solutions B^ = 3C and B^ = 12C, and the solution we are seeking 
is given by B = \/^C 1.73VC. Thus, we find that the curve e = e{p) encounters a 
singularity at B = \/3C and eq. (A. 4) definitely has no solutions for B > VSC. In fact, 
even before reaching this singularity there is a point where the curve becomes multiple- 
valued. At this point the tangent to the curve dpS goes to infinity. Thus by considering 
the system 

p = 0, {dper^ = 0, (92e)"' = 0, (A.24) 
we find using similar methods as above that 

45^2-72Si°C-565^C^-h3220S^C^-28627S''C*^-M09308S^C^-138384C^ = 0. (A.25) 

The solution of interest for us is — \/V3\JlC ~ 1.67\/C', thus there are no solutions 

to eq. (A. 4) for B > -Bmax = \/5 — \/T3y/2C. This situation is depicted in Figure 3. 

Summary: We end the appendix by summarizing our findings. To find solutions with 
the right boundary conditions then C > 0. For B < and < B < -Bmax, where 
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-Smax = \/5 — \/l3\/2C, the sohitions to the quintic equation are well approximated by 
the solutions to the cubic equation with at most 0(10"®) corrections. At i? ~ -Bmax the 
solution to the cubic approximates the solution to the quintic with corrections as big as 
up to 40%. Beyond B > Smax there are no solutions to the quintic equation satisfying the 
specified boundary conditions. 
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